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Abstract 

We study the thermal properties of a pinned disordered harmonic chain weakly perturbed by 
a noise and an anharmonic potential. The noise is controlled by a parameter A — s- 0, and the 
anharmonicity by a parameter A' < A. Let k be the conductivity of the chain, defined through 
the Green-Kubo formula. Under suitable hypotheses, we show that k = 0{X) and, in the absence 
of anharmonic potential, that k ~ A. This is in sharp contrast with the ordered chain for which 
K ^ 1/A, and so shows the persitence of localization effects for a non-integrable dynamics. 
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1 Introduction 



The mathematically rigorous derivation of macroscopic thermal properties of solids, starting from their 
microscopic description, is a serious challenge [5J[T7|. On the one hand, numerous experiments and 
numerical simulations show that, for a wide variety of materials, the heat flux J is related to the 
gradient of temperature T through a simple relation known as Fourier's law: 

J = -k(T)VT, 

where k[T) is the thermal conductivity of the solid. On the other hand, the mathematical understand- 
ing of this phenomenological law from the point of view of statistical mechanics is still lacking. 

A one-dimensional solid can be modeled by a chain of oscillators, each of them being possibly 
pinned by an external potential, and interacting through a nearest neighbors coupling. The case of 
homogeneous harmonic interactions can be readily analyzed, but it has been realized that this very 
idealized solid behaves like a perfect conductor, and so violates Fourier's law [H]. To take account of 
physical observations, one thus needs to consider more elaborate models, where ballistic transport of 
energy is broken. Here are two possible directions. 

On the one hand, adding some anharmonic interactions can drastically affect the conductivity of 
the chain [31 [T2] . Unfortunately, the rigorous study of anharmonic chains is in general out of reach, and 
even numerical simulations do not lead to completely unambiguous conclusions. In order to draw some 
clear picture, anharmonic interactions are mimicked in [21 E] by a stochastic noise that preserves total 
energy and possibly total momentum. The thermal behavior of anharmonic solids is, at a qualitative 
level, correctly reproduced by this partially stochastic model. By instance, the conductivity of the one- 
dimensional chain is shown to be positive and finite if the chain is pinned, and to diverge if momentum 
is conserved. 

On the other hand, impurities constitute another element that can affect the conductivity of an 
harmonic chain. In |22| and |10[ , an impure solid is modeled by a disordered harmonic chain, where the 
masses of the atoms are random. In these models, localization of eigenmodes induce a dramatic fall off 
of the conductivity. In the presence of everywhere onsite pinning, it is known that the chain behaves 
like a perfect insulator (see Remark 1 after Theorem [T]). The case of unpinned chain is more delicate, 
and turns out to depend on the boundary conditions [Hj. The principal cases have been rigorously 
analyzed in [24] and [1]. 

The thermal conductivity of an harmonic chain perturbed by both disorder and anharmonic inter- 
actions is a topic of both practical and mathematical interest. We will in the sequel only consider a 
one-dimensional disordered chain with everywhere onesite pinning. Doing so we avoid the pathological 
behavior of unpinned one-dimensional chains, and we focus on a case where the distinction between 
ordered and disordered harmonic chain is the sharpest. We will consider the joint action of a noise 
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and an anharmonic potential ; we call A the parameter controlling the noise, and A' the parameter 
controlling the anharmonicity (see Subsection 12.11 below) . 

The disordered harmonic chain is an integrable system where localization of the eigenmodes can be 
studied rigorously i6, • However, if some anharmonic potential is added, very few appears to be known 
about the persistence of localization effects. In |13) . it is shown trough numerical simulations that an 
even small amount of anharmonicity leads to a normal conductivity, destroying thus the localization of 
energy. In |20j , an analogous situation is studied and similar conclusions are reached. This is confirmed 
rigorously in 'S', if the anharmonic interactions are replaced by a stochastic noise preserving energy. 
Nothing however is said there about the conductivity as A — !• 0. Later, this partially stochastic system 
has been studied in |12j . where numerical simulations indicate that k ^ A as A — > 0. 

Let us mention that, although the literature on the destruction of localized states seems relatively 
sparse in the context of thermal transport, much more is to find in that of Anderson's localization and 
disordered quantum systems (see [3] and references in [31II2])- There as well however, few analytical 
results seem to be available. Moreover, the interpretation of results from these fields to the thermal 
conductivity of solids is delicate, in part because many studies deal with systems at zero temperature: 
the time evolution of an initially localized wave packet. 

The main goal of this article is to establish that disorder strongly influences the thermal conductivity 
of a harmonic chain, when both a small noise and small anharmonic interactions are added. We will 
always assume that A' < A, meaning that the noise is the dominant perturbative effect. Our main 
results, stated in Theorems [T] and [5] below, are that k ~ 0{X) as A ^ 0, and that k ^ A if A' = 0. 
Stricly speaking, our results do not imply anything about the case where A' > and A = 0. However, in 
the regime we are dealing with, the noise is expected to produce interactions between localized modes, 
and so to increase the conductivity. We thus conjecture that k = 0{X') in this later case. This is in 
agrement with numerical results in [2Q^, where it is suggested that k coud even decay as e^^/^ for 
some c > 0. 

In the next chapter, we define the model studied in this paper, we state our results and we give 
some heuristic indications. The rest of the paper is then devoted to the proof of Theorems [T] and [H in 
a way described after the statment of Theorem [5] 

2 Model and results 
2.1 Model 

We consider a one-dimensional chain of N oscillators, so that a state of the system is characterized by 
a point 

X = {q,p) = (qi,...,gAr,Pi,...,pjv) e M^^, 
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where qk represents the position of particle k, and pk its momentum. The dynamics is made of a 
hamiltonian part perturbed by a stochastic noise. 



The Hamiltonian. The Hamiltonian writes 
H{q,p) = Hhar{q,p) + X'Hanh{q,p) 

^ N N 

= ^^{pl + ''kql + {Qk+i - qkf) + A' ^ (t/ {qk) + V{qk+i - qk)) , 
k=l fe=l 

with the following definitions. 

• The pinning parameters Uk are i.i.d. random variables, which law is independent of N. One 
assumes this law to admit a bounded density, and that there exist constants < v- < 1/+ < oo 
such that 

P(l/_ < Uk < u+) = 1. 

• The value of gjv+i depends on the boundary conditions (BC). For fixed BC, one puts qN+i = 0. 
For periodic BC, one puts gjv+i = qi- For further use, one also defines qo = qi for fixed BC, and 
qo = 9jv for periodic BC. 

• One assumes that A' > 0. The potentials U and V are symmetric, meaning that U{—x) = U{x) 
and V{—x) = V{x) for every a; e M. One assumes that U and V belong to C^^p(M), the space 
of infinitely differentiable functions with polynomial growth. One assumes moreover that 



/ 

Jr 



e-^^^Ma; < +00 and d^U{x) > 0, 
and that there exists c > such that 

c < l + X'dlV{x) < c-^. 

For x,y G M**, let {x,y) = xiyi + • • • + XdVd be the canonical scalar product of x and y. The 
harmonic hamiltonian Hhar can also be written as 

Hhar{q,p) = ^{p,p) + ^{q,^q), 

if one introduces the symmetric matrix $ S R^^^ of the form $ = —A + W, where A is the discrete 
Laplacian, and W a random "potential". The precise definition of $ depends on the BC: 

$j,fe = (2 + Uk)5j^k - Sj,k+i - Sj,k-i (fixed BC), 

^j,k = + J^k)Sj,k - Sj,k+i - Sj^k-i - Si,N - Sn,i (periodic BC), 

for l<j,k<N. 
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The dynamics. The generator of the hamiltonian part of the dynamics is written as 




with 



N 



and 




{dxU{qk) + da;V{qk - qu-i) - d^^Viqk+i - qk)) ■ dj 



The generator of the noise is defined to be 



N 



XSu = X^{u{...,-pk,..-)-u{...,pk,...)) 



fe=i 



with A > A'. The generator of the full dynamics is given by 



L = A + XS. 



We denote by X*^ symply by X*-{x), the value of the Markov process generated by L at time 

t>0, starting from x = {q,p) G M?^ . 

Expectations. Three different expectations will be considered. We define 

• /xt: the expectation with respect to the Gibbs measure at temperature T, 

• E: the expectation with respect to the realizations of the noise, 

• E^: the expectation with respect to the realizations of the pinnings. 

We will also write Et for E^t- 

The Gibbs measure /xt is explicitly given by 



where Zt is a normalizing factor such that fiT is a probability measure on M?'^ . We will need some 
properties of this measure. Let us write 




^-ie-^(-)/r 



p'{pi)...p'iPN).p"{q) 



with p'{pk) = e-p'^/2'^/\/2^ for 1 < fc < AT. 



When A' = 0, the density p' is Gaussian: 



p'{q) = (27rT)-^/2.(det$)i/2.e-<«'*«>/2T_ 
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Since Vk > i'- > 0, it follows from Lemma 1.1 in [5| that |($ < Ce '^l-' for some constants 

C < +00 and c > independent of N. This implies in particular the decay of correlations 

When A' > 0, the density p' is not Gaussian anymore. We here impose the extra assumption 
that is large enough. In that case, our hypotheses ensure that the conclusions of Theorem 3.1 in 
[7] hold: there exist constants C < +00 and c > such that, for every /, g G C^„^p{R^) satisfying 
Mt(/) = ^J■Ti9) = 0, one has 

\Mf-9)\ < Ce-='^(S(^)-S(«))(MT((V,/,V,/)).A.T((V,g,V,.g)))'^'. (2.1) 

Here, S(u) is the support of a function u, defined as the smallest set of integers such that u can be 
written as a function of the variables xi for / G S(u), whereas d(S{f),S{g)) is the smallest distance 
between any integer in S(/) and any integer in S(.g). Using that /iT('Zfe) = for 1 < k < N ^ one 
checks from (|2.1[) that every function u G C~„p(M^) with given support independent of N is such that 
||u||Li(^y) is bounded uniformly in N. 

The current. The local energy of atom k is defined as 

Gfc — ^k,har ~1~ A C^^anh 

with 

ek,har = Y + + - (lk-1? + \{qk+i - qk? for 2 < fc < iV - 1, 

and 

TTI \ I - qk-l) , V{qk+1 - qk) f O ^ U ^ AT 1 

ek,anh = U(qk) H \ for 2 < k < N - 1. 

For periodic B.C., these expressions are still valid when fc = 1 or fc = A^. For fixed B.C. instead, one 
multiplies by a factor 2 all the terms involving the differences (go — qi) or [qN+i — qN) in the previous 
expressions. These definitions ensure that the total energy is the sum of the local energies. 
The definition of the dynamics implies that 

defe = {jk-i - jk) dt 

for local currents 

jk — jk,har ~t~ A jk,anh 

defined as follows for < fc < iV. First, for 1 < fc < TV - 1, 

jk.har = 7:iPk + Pk+l){qk - qk+l) and jk.anh = T^iPk + Pk+l) O^V {qk - qk+l) ■ (2.2) 
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Next, jo.i = jN,N+i = for fixed B.C. Finally, jo and are still given by (|2.2I) for periodic B.C., 
with the conventions po — pn and pn+i — Pi- The total current and the rescaled total current are 
then defined by 

N N 

Jn = JNJiar + X' JN,anh = ^ Jfc.h.ar + A' ^ jfc^an/i (2.3) 

k=l k=l 

rr i \ ' T J N .har ^/JN.anh In A\ 

Jn — JN,har + ^ JN,anh — ' 



2.2 Results 

For a given realization of the pinnings, the (Green-Kubo) conductivity k = k(A, A') of the chain is 
defined as 

k(A,A') = 7^1™ 1™ nt,N[\\') = lim lim Et \ ^ f JNoXf^y.ds] (2.5) 

if this limit exists. One expects the choice of the boundary conditions to play no role in this formula 
since the volume N is sent to infinity for fixed time. The homogenized conductivity is defined by 
replacing by Ej^Et in (j2.5p . By ergodicity, one expects that the conductivity and the homogenized 
conductivity coincide for almost all realization of the pinnings (see [5 ). The dependence of k(A, A') 
on the temperature T will not be analyzed in this work, so that one can consider T as a fixed given 
parameter. 

We first obtain an upper bound on the homogenized conductivity. 



Theorem 1. Let < A' < A. With the assumptions introduced up to here, if is large enough, and 
for fixed boundary conditions, 

2 

as A ^ 0. (2.6) 



1 ( I r Y 

— hmsuphmsupEi^ET^ / JNoXf^y.ds] = 0{X) 
t^oo N^oc \Vt Jo J 



Remarks. 1. When A = 0, the proof (see Section actualy shows that 

^limsupEi^^T^-^y JAroX(Q o)ds^ = ^(^"^) as t ^ oo. 

This bound had apparently never been published before. It says that the unperturbed chain behaves 
like a perfect insulator: the current integrated over arbitrarly long times remains bounded in L?{Ei,fj,T)- 

2. The proof (see Section [3]) shares some common features with a method used in [18] to obtain a weak 
coupling limit for noisy hamiltonian systems. In our case, one can indeed see the eigenmodes of the 
unperturbed system as weakly coupled by the noise and the anharmonic potentials. 

3. The choice of fixed boundary conditions just turns out to be convenient for technical reasons (see 
Section HI). 



7 



4. The hypothesis that is large enough is only used to ensure the exponential decay of correlations 
of the Gibbs measure when A' > 0. 

Next, in the absence of anharmonicity (A' = 0), one gets more refined results. 

Theorem 2. Let A > 0, let A' = 0, and let us assume that hypotheses introduced up to here hold. For 
almost all realizations of the pinnings, the Green-Kubo conductivity (|2.5p of the chain is well defined, 
and in fact 

0) - ^ /iS. A- E^^- (t! I ° ^(>^o) ' (2.7) 

this last limit being independent of the choice of boundary conditions (fixed or periodic). Moreover, 
there exists a constant c > such that, for every A g]0, 1[, 

cA < k(A,0) < c-^X. (2.8) 

This article is devoted to the proof of these theorems, which is constructed as follows. 
Proof of Theorems [7] and The upper bound (|2.6I) is derived in Section [3j assuming that Lemma [1] 
holds. This lemma is stated and shown in Section!?]; it encapsulates the informations we need about 
the localization of the eigenmodes of the unperturbed system (A = A' = 0). The existence of k(A, 0) for 
almost every realization of the pinnings, together with (j2.7p . are shown in Section [5] Finally, a lower 
bomid on the conductivity when A' = is obtained in Section [S] This shows (|2.8p . □ 

2.3 Heuristic comments 

We would like to give here some intuition on the conductivity of disordered harmonic chains perturbed 
by a weak noise only, so with A > small and A' — 0. We will develop in a more probabilistic way 
some ideas from p^2^. Our results cover the case where the pinning parameters Vk are bounded from 
below by a positive constant, but it could be obviously desirable to understand the unpinned chain as 
well, in which case randomness has to be puted on the value of the masses. We handle here both cases. 

Let us first assume that Vk > c for some c > 0, and let us consider a typical realization of the 
pinnings. In the absence of noise (A — 0), the dynamics of the chain is actually equivalent to that of N 
independent one-dimensional harmonic oscillator, called eignmodes (see Subsection 14.11 and formulas 
(j4.5M.6D in particular). Since the chain is pinned at each site, the eigenfrequencies of these modes are 
uniformly bounded away from zero. As a result one expects all the modes to be exponentially localized. 
One can thus naively think that, to each particle, is associated a mode localized near the equilibrium 
position of this particle. 

When the noise is turned on (A > 0), energy starts being exchanged between near modes. Assume 
that, initially, energy is distributed uniformly between all the modes, except around the origin, where 
some more energy is added. We expect this extra amount of energy to diffuse with time, with a variance 
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proportional to k(0, A) . i at time t. Since flips of velocity occur at random times and with rate A, one 
could compare the location of this extra energy at time t to the position of a standard random walk 
after n = Xt steps. Therefore, denoting by Sk the increments of this walk, one gets 

1 " 2 
K{X,0) ^ {{^T5,) ) ^ A. 

^ k=l 

This intuitive picture will only be partially justified, as explained in the Remark after the proof of 
Theorem [T] in Section [31 

Let us now consider the unpinned chain. So we put ui^ — and we change by p^/ruk in the 
Hamiltonian, where are i.i.d. positive random variables. We consider a typical realization of the 
masses. In contrast with the pinned chain, the eignefrequencies of the modes are now distributed in 
an interval of the form [0, c], for some c > 0. This has an important consequence on the localization of 
the modes. One expects indeed that the localization length I of a mode and its eignefrequency uj are 
related through the formula I ^ XjuP'. 

Here again, the noise induces exchange of energy between modes, and one still would like to com- 
pare k(A,0) .< with the variance of a centred random walk with increments 5k- However, due to the 
unlocalized low modes, bk can now take larger values than in the pinned case. Assuming that the 
eigenfrequencies are uniformly distributed in [0, c], one guesses that, for large a, 

P(|4|>a) - P(\lu?>a) ^ 

This however neglects a fact. Since energy does not travel faster than ballistically, and since successive 
flips of the velocity are spaced by time intervals of order 1/A, it is resonable to introduce the cutt-off 
P(|(5fc| > 1/A) = 0. With this distribution for |(5fe|, and with n — Xt, we now find 

^ k=l 

This scaling is numerically observed in [12] . The arguments leading to this conclusion are very approx- 
imative however, and it should be desirable to anlyze this case rigorously as well. 

3 Upper bound on the conductivity 

We here proceed to the proof of Theorcm[T] We assume that Lemma[T]in Section|3]holds: there exists a 
sequence {un)n>i C L^(E^/xt) such that -AharUN = JN.har, and that iuN)N>i and {AanhUN)N>i are 
both bounded sequences in L^(Ei,/iT)- MoreoverwAr is of the form ujv(<?,p) — {q,aNq) + {p,jnP) +cn, 
where a,7 G are symmetric matrices, and where cn G K. 

Proof of (|2.6p . Let < A' < A, and let un be the sequence obtained by Lemma[T]in Section^ Before 
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starting, let us observe that, due to the special form of the function u^, one has 

N 

AanhUN = ^0/(9) Pi, (3.1) 

1=1 

with 

AT 

(j)i{q) = 2^jk,i ■ {dxV{qk+i - Qk) - dxV{qk - Qk-i) - dxU{qk)j, (3.2) 
fc=i 

where (7fe,i)i<fe,/<Af are the entries of 77V. It follows in particular that 

AanhUN = -^{-S)AanhUN- (3.3) 

Now, since J'n = JN,har + ^' JN,anh, ouc has, by Cauchy-Schwarz inequality, 

E,Et ^ ° =^ E.Et (^-^ ^ jNMar « +(A')'E,Et (^"^ ^ JiV^anh O X, 

+ 2A' I^E.Et (^-^ ^ Jiv,/.ar o ds^ . E.E/iT ^ J}v,a„/. o X' ds^ j . 
Since JN.anh = \{—S)JN.anhi a classical bound 15^ furnishes 

Ej/Et y Jn ,anh ° ds^ < C EvI^t{Jn ,anh ■ { — ^S) ^jN,anh) < '^^v f^T {Jn ^anh) 

where C < +00 is a universal constant. By (12. ip . ^vIJ-t^Jn anh) bounded uniformly in A^. Therefore 
limsuplimsupE^^ET ( ^ / JN,anh ° X'' ds\ = C'(A"^). 

t-i-oo N=foa \\t Jo J 

It suffices thus to establish that 

limsuplimsup Ei^E^T ( ^ / >JN,har o X" dsj = 0{X). 

t^oo Ar->oo Wt Jo ' J 

One writes 

/ A' \ 

JNJiar = -AharUN = -Xwat + A'^a„;iMAr + AS'tiAT = -iUAT + A S'^/d - — Ann,, j MAT, 

where the second equality is obtained by means of (j3.3p . Therefore 

y Jnmt °X'ds ^ ^ J LuN o X" ds + ^ J S (^Id- ^Aanh^ "jv ° X' ds 

1 _ uoX^-u A 
= —Mt p 

where A^t is a martingale given by 



^ + ^J^s(^Id- ^Aan,^ UN O ds, (3.4) 



Mt = f ^S'jUwoX, (dNi-Ads), 
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with the Poisson process that flips the momentum of particle j. 

It now suffices to establish that the three terms in the right hand side of <\'SA\i are 0{\) in L'^{Ei,Et). 
Up to the end of this proof, we use the simpler notation (•), for ^t{')- Let us first show that 
{un.{~S)un) < 4||ujv||l2(^^)- Writing 

UN = Wj^ + W}v + ""at + CAT 



with 



one finds indeed 



u 



Pj,.Q ^ ^ -yijPiPj, u!"/'^ = ^ lt,iPh u%'' = {q, ONq), 



{uN ■ {-S)un) = {uN ■ = 4 . -u^^'") 

and 

{UN . um) = {uT' • uT") + + + cn)') + 2 {{uP/'' + ujf + cn) . uPf'"). 

The claim follows since ((w^^'"^ + m^'' + cn) ■ w^^'") = 0. 
So first, 



1 



|2 



Next, 

Finally, by a classical bound [TS] , 



V 2 



Et ^ 5 (^W - ^^a„h) WW O X, ds^ <CX^(^S (^Id - ^^a„h) UN, {-XS)-^S (^Id - ^^a„h) UN 

= C\(^ ^/d - ^Aa„/i^ UN, (-5) - ^Aa„/i j UJV^ 

= CA ^(WAT, (-5)wAr) + ^(y)^Pa„/,Ujv||L2(^^)^ 
< CA (^||wAr||^2(^^) + Pa«ftUAr||L2(pr)) 

where one has used (|3.3p and 

JV 

(uat, Aan/iUAf) = (((g, ajyg) + {p,jNP) + Cn) ■ y^0;(g)p;^ = 

to get the second equality. Taking the expectation over the pinnings, the proof is completed since 
{un)n and {AanhUN)N are bounded sequences in L^{E^^t)- D 
Remark. When A' — 0, formula p.4p becomes 

^ fjNoX^ds = -1m- "^°^''"^ +4= fsuNoX^ds. (3.5) 
\/iyo \/i Vi Vi Jo 
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Now, since Sun = -4:J2i<k^i<N IkjPkPi, one has 

J'jyoXt-sds = — / Jjv o ds and / SujyoXt-sds = / SuN°Xsds. 
Jo Jo 

The measure on the paths being invariant under time reversal, one has 

r \ 

Jjq o Xgds . \ Sun o Xgds \ =0. 



Therefore, from (|3.5p . one deduces that 



N 2 / \ 2 

t \ ,1.9 / \ rt \ 



where r{t) is quantity that vanishes in the limit t ^ oo. We see thus that our proof does not completly 
justify the heuristic developed in Subsection 12.31 due to the second term in the right hand side of this 
last equation. As explained after the statment of Lemma [1] below, the sequence ujv should not be 
unique. It could be that a good choice of sequence un should make this second term to be 0{X^). 



4 Poisson equation for the unperturbed dynamics 

In this section, we state and prove the folowing lemma. Fixed BC are assumed for the whole section. 

Lemma 1. Let X' > 0, and assume fixed boundary conditions. For every N > 1, and for almost every 
realization of the pinnings, there exist a function un of the form 

UN{q,p) = {q,aNq) + {p,jNP) + CN, 
where unj'Jn G B.^^^ are symmetric matrices and where cm S such that 

— AharUN = JN,har- (4.1) 

Moreover, the functions u m can he taken so that 

{un)n>i and {AanhUN)N>i are boundcd scquenccs in L^{E^fiT)- 

Remarks. 1. The parameter A' only plays a role through the definition of the measure ^t- 
2. For a given value of iV and for almost every realization of the pinnings, the unperturbed dynamics 
is integrable, meaning here that it can be decomposed in N ergodic components, each of them corre- 
sponding to the motion of a single one dimensional harmonic oscillator (see Subsection 14.11 and (j4.5M.6p 
in particular). This has two implications. First, since (j4.ip admits a solution, we conclude that the 
current JTat is of mean zero with respect to the microcanonical measures of each ergodic component 
of the dynamics. Next, the solution un is not unique since every function / constant on the ergodic 
components of the dynamics satisfies —A^arf = 0. 

Proof of Lemma[l\ To simplify notations, we will generally not write the dependence on N explicitly. 
The proof is made of several steps. 
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4.1 Identifying (ttAr)iv>i: eigenmode expansion 

Let z > and let 1 < ^, m < A^. Let us consider the equation 

(z - Aanh)ui,m,z = qiPm- 

The solution ui^m,z exists and is unique. It is given by 

/"OO 

ui,mA^) = / e"''*.q/oX?o(,N(x).p„ioXfoo)(a;)ds. (4.2) 
Jo 

We will analyze ui,,n,z to obtain the sequence un- Although we assumed fixed BC, all the results of 
this subsection apply for periodic BC as well. 

Solutions to Hamilton's equations. The matrix $ is a real symmetric positive definite matrix in 
M^^^, and there exist thus an orhonormal basis {£,'')i<k<N of R^, and a sequence of positive real 
numbers (a;^)i<fc<Ar, such that 

One can show that 

min{i^j : 1 < j < A^} < col < maxji/j : I < j < N} + A (4.3) 

for 1 < fc < A^. According to Proposition ILl in [16], for almost all realization of the pinnings, none of 
the eigenvalue is degenerate: 

Uj^ujk if j^fc, l<j,fc<iV. (4.4) 

In the sequel, we will assume that (|4.4p holds. 
When A = A' = 0, Hamilton's equations write 

dq — pdt, dp — —^qdt. 

For initial conditions {q,p), the solutions write 

N 

k=i ^ '^'^ ^ 

N 

^^{-^k{q,£.'')smujkt+{p,^'')cosuJkt)i''. (4.6) 

k=l 

An expression for ui,m,z- To determine ui^m.z, one just needs to insert the solutions (j4.5M.6p in 
the definition (j4.2p , and then compute the integral, which is a sum of Laplace transforms of sines and 
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cosines: 

N 



fe=i -"fe 

,2 



+ 5^ (Z,eKm,C'=)(-^(g,e)(9,^'=) + -^(p,^^)(p,^'=)) 

l<j^k<N ^^J-'^k ^j-^k ' 

+ 0{z\ (4.7) 

where (j, ^*^) denotes the f"^ component of the vector ^ , and where the rest term 0[z) is a polynomial 
of the form (g, ag) + (g, ^p) + (p, 7p) , where a and 7 can be taken to be symmetric. One defines 

2->-0 

One observes that ui^m is of the form (g, ag) + (p, 7p) where a and 7 can be taken to be symmetric. 
Defining the solution uj^. For fixed BC, the total current is given by 

JV-l 



2 

fe=i 



Setting 



for 2 < Z < iV and 



we define 



Wl = Un,N - - ht{un,n - 



JV 

One observes that un is of the form un = {q,aNq) + {PjJnP) + cn, where ajv and 7jv are symmetric 
matrices, and whee cjv G M. 

Let us show that un solves —AanhUN = Jn- One may assume that cat = without loss of 
generality. The current J7}v can be written as J-^ = {q,Bp). The function un has been obtained 
as the limit as z ^ of the function of the form = (g, cc^g) + (g, PzP) + (p, 'JzP) which solves 
{z — Ahar)uz = -Jn, and with and 7^ symmetric. Since 

{z - Ahar)uz = (g, {zaz + I3z^)q) + {q, [z^z - '^■{oiz - ^lz))v) + (P, (-2 - ^z)v). 

one has 

zoLz + \{{iz^ + */3l) = 0, zfiz - 2{,az - ^iz) = B, z-\{pz+ f^l) = 0. 

One knows that {az,Pz,lz) (Q!jv,0,7jv) as 2: 0, with ajv and7jv symmetric, so that — 2(a — $7) = 
B. Taking into account that ajv and 7jv are symmetric, one gets 

$7jv-7iv$ = ^{B-B^), 2aN = 2<i>jN - B. (4.8) 
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It is checked that, if two symmetric matrices and 7jv satisfy these relations, then un = {q,aNq) + 
{p,jNP) solves the equation —AharUN = Jn- 

4.2 A new expression for wi 

For 1 < Z < A^, let us write 

wi = {q,a{l)q) + {p,j{l)p) + c{l), 

where a{l) and are symmetric matrices, and where c{l) e K. A relation similar to (j4.8l) is satisfied: 
defining the matrices 

(Bil))^^^ = ki-iim,n)-Si^i.i{m,n) {2 < I < N) and {B{1)) = SN,N{m,n) ~ S,,i{m,n), 
for 1 < rrijU < N , one has 

2a{l) = 2$7(0-S(/), (4.9) 

for 1 < I < N. Therefore the knowledge of the matrices 7 impies that of the matrices a. 

An expression for the matrices 7(Z) can be recovered from (|4.7p with z — 0. We will now work this 
out in order to obtain a more tractable formula. We show here that, for 2 < Z < TV, 

N N 

j=i fc=i 

l-l N 

j=i k=i 

l-l N 

j=i k=i 

N N 

= -EEO'^')'(^'^')(^'^')' l<s^t<N (4.10) 
j=l fc=l 



and 



1 ^ I AT ckxl 1 /I <:k\2 

7m(1) = tE^^(^'^')(*'^')-iE^t(^'^')(^'^')' 1<^,^<^^- (4.11) 

k=l k=l *: 

Formula (j4.1ip is directly derived from (j4.7p . noting that 7(1) is the only symmetric matrix such 
that ?i;i(0,p) — ( 7s,t(l)psPt- To derive (j4.10p . one observes that 7(Z) is the only symmetric matrix 
such that wi{0,p) — ^ js.tiOPsPt- Stating from (|4.7p . one deduces 

— \ / UJ - — LO, 

l<j^k<N J « 

For fixed BC, the eigenvectors satisfy the following relations for 1 < j < N: 

(C^O) = {^^,N + 1) ^0 (by defintion), 
-(e,m-l) + (2 + i.„0(e,m)-(e,m + l) =u:^{i^,m), l<m<N. 
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So the following recurrence relation is satisfied: 

{e,m + l) ^ {2 + Vra-^]){e,m)-{i\m-l), l<m<N. (4.12) 
Let us first compute W2(0,p). Using (j4.12p . it conies 

(2, e') - (1, e)(2, e) ={2+^,- c.|)(i, ^^■) (i, e) - (2 + - c.^)(i, ^^■) (i, e) 

^ -{oj]~uji){i,i=){i,e)- 

Therefore 

AT 

l<j^k<N l<j,k<N k=l 

N 

= -(I,p)2 + E(l,a2(p,e')', (4.13) 
fe=i 

where the last equality follows from the fact that {^'')k forms an orthonormal basis. 
Let us now compute w/(0,p) for 2 < I < N . Again by (I4.12p . 

ii e") {I - 1, e) -{I- 1, e) a, e) = ((2 + lyi-i - - 1..^) -{1-2, e^)) {I - 1, e) 

- {i-h e) ((2 + i^i-i - Loi){i - 1, e'^) - - 2, e)) 
+ {i-i,^^){i-2,e)-{i-2,e){i-i,e)- 

Therefore 

wi{0,p) = - {l-l.e){l-hC''){p,e){p,^') + m-iiO,p) 

l<j^k<N 

N 

= - {I- i,pf + - i,e)'{p,ey + wi^,{o,p). (4.i4) 
fe=i 

Combining (j4.13p and (|4.14[) . one arrives to an expression valid for 2 < Z < iV: 



Wl 



3=1 \k=l ) 



Let us now write {j-.p)'^ = p^ and 
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One gets 



i-l N l-l 

'"1 



s,t ] = 1 k=l ] = ! 

l-l I l-l N 

s=l \j=l k=l 

N l-l N 
s=l j=l k=l 

l-l N 

+ E P.ftEEo-'^')'(-'^')(^'^')- 

l<s^t<N j = l k=l 

In this formula, the coefficients of coincide with 7s,s(0 given by (j4.10p for Z < s < A^, and the 
coefficients of PsPt with s t coincide with the first expression of 7s,t(0 given by (j4.10p . To recover 
the coefficients js.sil) for 1 < s < / — 1, just use the fact that (^'^)fc and {\k))k are orhonormal basis: 

l-l N N / N \ 

EEo-'^')'(-'^')'-i = E i-Eo-'^')M 

j = l fe = l y j = i J 

N N N N N 

= E(-i^')'-i-EE(j'^')'(-'^'^')' = -EEo-'^')'(-'^')'- 

fc=l j=l k=l j=l k=l 

The second expression for the coefficients js,t{l) with s 7^ t in (|4.10l) is obtained by a similar trick. 
4.3 Exponential bounds 

We show here that there exist constants C < +00 and c > independent of N such that 

E.(4,(0) < Cexp(-c(|j-Z| + |fc-;|)), E.(7|,fe(0) < C exp ( - c(|j - ;| + |fc - Z|)) , (4.15) 

for 2 < ^ < and for I < j,k < N. This is still valid for Z = 1 if one replcaces |A: — Z| by min{|fc — 
1|, j/c — A^l} and |j — Z| by min{|j — 1|, \j — A^l}. Due to (14. 9p . it suffices to establish these bounds for 
the matrices 7. 

Let us first observe that the almost sure bounds 

|7..(0I < 1 (2 < / < iV), |7.,(1)| < 

hold for 1 < s, i < A^. This is directly deduced from (|4.10l) and (|4.11l) by taking absolute values inside 
the sums if needed, using the bound < using that (^'=)fc and {\k))k are orhonormal basis, 

and Cauchy-Schwarz inequality if needed. By (j4.3p . minjo;^ : 1 < fc < N} > c > 0, where c does not 
depend on N. In particular Ej, (|7s.t|P) < CpE^|7s_j| for every p > 1, so that we only need to bound 

Ei/|7s,t|- 
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We now will apply localization results originaly derived by Kunz and Souillard [TB], but we follow 
the exposition by [11,. From (|4.10p and (|4.1ip . we see that we are looking for upper bound on the 
absolute value of sums of the type 

k 

and of the type 

Y,{r,er{s,e){t,e), Y.('-^e){s,e)'{t,e), Y.('-^e){s,e){t,e)\ r<s<t. 

bmee I {r, I < 1 for 1 < r, fc < A^, all of them can be bounded by 

k 

By the formula before Lemma 4.3 in [11], and the lines after the proof of this lemma, one concludes 
that there exist constants C < +oo and c > independent of N such that 



k 

Together with the remarks formulated up to here, this alows to deduce (|4.15p . 

4.4 Concluding the proof of Lemma [T] 

One writes 

E^firiujf) = ^EE^^t(w„.W;„) and Ei./iT ((^an/i.WAf)^) = E^^T(^a«feWm ■ AgnhWn)- 

m.n m,n 

One will establish that there exist constants C < +00 and c > such that 

\E,fiT{w„,.Wn)\ < Ce-^l"-"l and \E,^iT{AanhW„, . AanhWn)\ < Ce-=l™-"l (4.16) 

for 1 < m,n < N. This will conclude the proof. 

Let us fix 1 < m,n < N. Let us first consider |E^/iT(wm ■Wn)\- Let us observe that the functions 
wi are of zero mean by construction, and so the relation 

J2 '^j.j(^) / + J2 "j^fe(^) / ^^^^ + ^(^) = (4.17) 

j j,k 

holds for 1 < / < A^. Using this relation, one computes 

= E a.,,j{m)ak,i{n) J qiqjqkqi df^T + ^ ai,j{m)-fk,i{n) j qiqjPkPid^T 

i^j.k.l i,j.k.l 

"fi,jim)ak,iin) J PiPjqkQi dfj.T + ^ li,j{m)jk,i{n) J piPjPkPi dfir - c{m)c{n) 



i,j,k,l i,j,k,l 

Si + S2 + S3 + S4 — c{m)c{n). 
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Using (|4.17p and the fact that J dfiT = T, the sum 5*1 is rewritten as 

^1 = ^ a^^j{m)ak,iin) J [q^q^ - j qiqj dfixj [quqi ~ J qkQi d/ir) d/Ltr 

i^j.k.l 

+ X! aij("i)afe,/(w) J qiq-jdfiT J qkQid^ir 

i,j,k,l 

= X! o;tj{m)akj{n) J (^qtqj - J qiq-j dfixj (^qkqi - J qkQl d^r^ d^r 

i.j,k,l 

+ r^^7i,i(™)7jj('^) +Tc{m)^^jij{n) + Tc(n) ^ 7i^i(m) +c{m)c{n). 

i,j i i 

Then, stih using (|4.17p . one obtains 

5'2 + 5'3 = -Tc(m) ^7i^i(n) - rc(n) ^7,;^i(m) - 2r^^7i,i (m)7j J (n). 

Finally, the terms in the sum 6*4 are non zero only when 

i=j = k = l, i = j,k = l,i ^ k, i = k,j = l,i ^ j, i = l,j = k,i ^ j. 
Using that /p^d^T = 3(/p^ dfir)^ and that Jp^ dfiT — T, S4 is seen to be equal to 

54 = (^7,,,(m)7j-j (n) + 27,^(771)7,^(71)^ 

Therefore 

^iT{wm-Wn) = ^ aij{m)ak,i{n) J {l^lJ ~ J QiQj dfJ-T^ (qkQi - J qkQi dfj,T^ dfj,T 

i,j,k,l 

+ 2T2^ 7,^(777)7,^(71). 

hi 

Applying the decorrelation bound (12.11) and the exponential estimate (|4.15p . one obtains the result. 
Let us next consider lE^^Ti^anhWrn ■ AanhWn)\- From p.ip . one has 

^J■T{Aa„hWm ■ AanhWn) = ^ HT {(l>tiq, m) (f>s{q, n) PsPt) = ^ (^t (g, m) (f>t{q,n)). 

s,t t 

Now, from (15^ . one has (t)t{q,k) = 7s,t(fc)/0s(<z) for 1 < fc < A^. Here ps{q) = ps(qs_i, g^, g^+i) is 
a function of mean zero since the potentials U and V are symmetric. One writes 

i^T{AanhWm ■ AanhWn) = Is.tjm) -fs'A^) i^riPs ■ Ps')- 

t s.s' 

Applying the decorrelation bound (j2.ip and the exponential estimate (j4.15p . one obtains the result. □ 

5 Convergence results 

In this Section we show the convergence result (|2.7p . We assume thus A > and A' — 0. 
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We start with some definitions (see [H] for details) . The dynamics defined in Section [3] can also 
be defined for a set of particles indexed in Z instead of TLjsi. Points on the phase space are written 
X — {q,p), with q = {qk)kez and p — {pk)kez- Let us denote by C the generator of this infinite- 
dimensional dynamics, and by fir the corresponding Gibbs measure at temperature T. We extend 
the definition (j2.2p of local currents jk to all fc € Z {jk — jk,har since A' = 0). li u ~ u{x^v), with 
v = {h'k)kez a sequence of pinnings, and if fc e Z, we write Tku{x, i>) = uirkx, Tkv), where 

{Tkq)j = qk+j, {Tkp)j = Pk+j, (Tkiy)j = i^k+j- 

We will use the hermitian product {u,v)n — P-t{u - v*), where v* is the complex conjugate of v, 
and where N £^ ov N — oo for the infinite dimensional dynamics. We will always write the index N to 
avoid confusion with the canonical scalar product (•, •) on W^. We denote by ^ •, • ^ the inner-product 
defined, for local bounded functions u and v, by 

E^((u,Tfci;)oo - p.t{u)p.t{v)) 

k£Z 

and by % the corresponding Hilbert space, obtained by completion of the bounded local functions. 

We start with two lemmas. We have no reason to think that Lemma [2] still holds if an anharmonic 
potential is added, and this is the main reason why we here restrict ourselves to harmonic interactions. 

Lemma 2. There exists a constant C < +cx) such that, for any realization of the pinnings, for the 
finite dimensional dynamics with free or fixed B.C., or for the infinite dynamics, for any fc > 1 and 
for any ^ G Z^v (resp. fc e Z for the infinite dynamics), one has 

l|i'j/||L^(M.) < c^- (resp. < C^;. 

Proof. Let us consider the infinite dimensional dynamics ; other cases are similar. We can take I = 
without loss of generality. The function jo is of the form jo = {q, aq) + {q, fip) + {p, 7p), with a = 7 = 
and P defined by 

Pi.j = ^('5o,o(j,i)-'5(i4)(i,j)-|-Jo,i(j,j)-'5i,o(i,i))- 

Now, if u is any function of the type u = (q, aq) + (g, (ip) + {p, jp) , then Cu = (q, a'q) + {q, f3'p) + {p, j'p) 
with 

where 7 is such that (7)4, i = and — 7ij- for i 7^ j. Thus 

^''jo = {q,o^ik)q) + {q,f3{k)p) + {p,i{k)p), 

and there exists a constant C < +00 such that Q,j — whenever |i| > Cfc or |j| > Cfc and such that 
Ci,j ^ otherwise, with C one of the three matrices a(fc), P(k) or 7(fc)- The claim is obtained by 
expressing \\C'ji\\j^2(^^^) in terms of the matrices a(k), fi{k) and "i{k)- D 

Explicit representation for the matrix in Lemma 1.1. in [S] allows to deduce the following lemma. 
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Lemma 3. Let f and g be two polynomials of the type {q, aq) + (p, /3p) + (p, 773), and assume that there 
exists 71 e N such that aij — — ji^j — whenever \i\ > n or \j\ > n. Then there exists c > such 
that, for fixed or periodic B.C., 

\{f,9)N-{f,9)oc\ ^ Oie-^"") as iV ^ cx). 

Let then D = {z e C : 3fiz > 0}. For every z G D, let be the unique solution to the equation 

= jo. (5.1) 

We know from Theorem 1 in [sj^ . and from its proof, that 

lim <C Uz,jo 3> exists and is finite (5-2) 

2^-0 



and that 



limz<Cu2,Uz^= 0. (5-3) 



For z € D and iV > 3, let Uk^z,N be the unique solution to the equation 

(z - L)uk,z,N = jk- (5-4) 

so that 

Uz,N ■= — ^y^Mfc,2,JV solves {z - L)uz.N = Jn- (5.5) 

^ k 

Lemma 4. For fixed or free boundary conditions and for almost all realizations of the pinnings, 
lim lim ^.t{uz,n -Jn) = lim lim E.^^it(uz.n ■ Jn) = lim < Wz, jo >, 
lim lim z ut(uz n -Uz n) — lim lim zEi,riT(uz n -Uz n) — 0. 

Proof. By (|5.2p and (|5.3p . it suffices to establish separately that, for every z e B, and for almost every 
realization of the pinnings, 

lim ijlt{uz,n -Jn) ^ <.Uz,3o'>, lim E^^itUz.n ■ Jn) =<W2,io>, 

lim z [iriuzN .Uzn\ = z<u^,u^>, lim zEv[1t{uz,n ■Uzn\ = z<$:uz,Uz:^. 

The proof of these four relations is in fact very similar, and we will focus on the first one. We proceed 
in two steps: we first show the result for |z| large enough, and then extend it to all z e B. 

First step. Here we fix z G B with |z| large enough. We first assume periodic boundary conditions. 
The function Uz solving (j5.ll) may be given by 



k>Q 



^ The model studied there is not exactly the same. The proof of the properties we mention here can be however 
readily adapted. 
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this series converging in virtu of Lemma [3] for \z\ large enough. Let now n > 1. One computes 

n oo 
fc=0 iGZ fc=n+l (GZ 

For every given fc, the sum over I is actually a sum over C k non-zero terms only, for some C < +00. 
From this fact and from Lemma [2j it is concluded that the second sum in the right hand side of (|5.6p 
converges to as n — >■ 00. Similarly one writes 



{Jn,u..,n)n = ^EE^"^''''^^-^'^*)^ + ]^EE^"^'^'^0-^'j*)^'- (5.7) 

s,t k=0 s,t k=0 

Here as well, the second term in (|5.7I) is such that 

n 

lim Iimsup-EE^^^'^'^(J-^'^'*)^ = 0- 
To handle the first term in (|5.7p , let us write 

n 

^"('^) = E E-^^^'^'^jO'^'j*)^- 

teZw fc=o 

One has in fact 

s,t /c=0 s 

The result is obtained by letting N 00, invoking Lemma [3] and the ergodic theorem, and then letting 
n 00. If one had started with fixed boundary conditions, then, for every fixed n, all the previous 
formulas remain valid up to some border terms that vanish in the limit — > 00 due to the factor 1/iV. 
Second step. Denote by £,n^^{z) and £(z) the complex functions defined on D by 

S.N.i'iz) = {uz^n,Jn)n and 2{z) =<u^,jo>. 

The first observation is that these functions are well defined and analytic on D. Moreover, similarly 
to what is proved in [5j, they are uniformly bounded on D by a constant independent of N and the 
realization of the pinning v. 

Let us fix a realization of the pinnings. The family {£,n,i' , > 1} is a normal family and by 
Montel's Theorem we can extract a subsequence {£'Nk,v}k>i such that it converges (uniformly on every 
compact set of D) to an analytic function /*. 

By the first step we know that f*{z) — £(2) for any real z > zq. Thus, since the functions involved 
are analytic, /* coincides with £ on D. It follows that the sequence {£jv,i/(-z)}Af>i converges for any 
z e © to £(z). □ 

Following a classical argument, we can now proceed to the 
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Proof of (|2.7I) . For any z > 0, one has 



1 /•< 1 /-t i-t 



^ I jNO-'^''ds = ^ / Luz,N ° X" ds + / Uz.AToX'^ds 

V Wo V Wo yt Jo 

= ^^-"'* — - + viL 

Here A^z,iv,f is a stationary martingale which variance is given by 

where one has used that ij-t{uz,n ■ ^harUz^w) — 0. Next one has 

' Uz.N ° X'' - Uz.n\'^ ^ 2 

and 



Et ■ —j < J ^^T{Uz,N -Uz^n) 



Et (^-^ J Uz,NoX''ds^ < z'^t ^it{uz.n -Uz.n)- 

Remimding that fiT{uz.N ■ {z — L)uz,n) — fJ'T{uz,N ■ Jn), the proof is completed by taking z — 1/t 
and invoking Lemma |31 D 

6 Lower bound in the absence of anharmonicity 

We here establish the lower bound in (j2.8p . and so we assume A > and A' — 0. We also assume 
periodic boundary conditions. We use the same method as in [S] (see also [12]). According to Section 
[51 it is enough to establish that there exists a constant c > such that, for almost every realization of 
the pinnings, for every z > Q and for every > 3, one has 

lJiT{jN-{z-LfjN) > C. (6.1) 

Proof of (|6.ip . For periodic B.C., the total current Jn is given by 

-^JV = 2 X! i^kPk+i - Qk+lPk) 

To get a lower bound on the conductivity, we use the following variational formula 

liT {Jn (z - L)-' Jn) = sup {2^t(Jjv . f) - /it if ■ (z - XS)f) - /xt (A^/ ■ (z - XS)-^ Aharf)} 

(6.2) 

where the supremum is carried over the test functions / £ C^^p{R^^). See [33] for a proof. We take 
/ in the form 

/ = a{q, Pp) with a G M and /3 = $Af, 

where M is the antisymmetric matrix such that Mij — Si_j-i — with the convention of periodic 

B.C.: Si^N+i = ^1,1 and 6n,o = Sn,n- 
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First, we have 

Aharf = a{p, I3p) - a{q, I3<^q) = a^P.^^p.pj 

since /3$ — $M$ is antisymmetric, and since = for 1 < j < A^. Since S{piPj) = —-ipiPj for i ^ j, 
one obtains 

{Af . (z - \S)-' Af) = -l^^^iA,,,,/ . A,,arf) = (l^lj + /^^J/^J-O 

for some constant C < +00. Next, since Spk — — 2pfe for < fc < A'', there exists some constant C < +00 
such that 



/xt(/ . (z - XS)f ) = a^T^z + 2A) ^ A,^ = a2T2(z + 2A)Tr [/3t$-i/3] 

< Ca2T2(z + 2A)A^. (6.4) 
Let us finally estimate the term ^t{ Jn /)■ 

fj-riJN f) ^ ^(^^Pi.jqiPj ■ Y {qkPk+i - qk+iPki^ ^ ^^^{p.^k+iiqtqk) ~ Pi,k{qiqk+i)) 

i,j A;GZjv i-,k 

= — E {{fi^^-%+i.k-{[i^^-')k+i,k) = — E {Mk,k+i-Mu+iM) 

= aT'^N. (6.5) 



By (|6.3p . (|6.4p . (j6.5p and the variational formula (|6.2p . we find that there exists a constant C < +00, 
independent of the realization of the disorder, of A and of N, such that for any positive a. 



1 



rMT(Jw {z - L)-^ Jn) > a- Ca^ (z + 2A) + 
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By optimizing over a, this implies 

^,MJm{z-L)-'Jn) > ^ (^(- + 2A) + 

This concludes the proof. □ 
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